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1. Introduction 

The current standard cosmological paradigm, the so called ACDM universe, establishes that 
the universe has a spatially flat geometry whose dynamics is nowadays dominated by dark energy 
with a minor contribution from the dark matter, both of them are of unknown nature and where 
neither radiation nor baryonic matter play an important role in the present dynamical state of the 
universe. In this scenario, structures grow as gravitational instabilities in the dark matter density 
field where the first objects are the smaller ones and subsequent mergers and accretion of the 
small objects on to the big ones lead to a hierarchical scenario of growth of structures. Galaxies 
are supposed to form as gas accretes and cools down in the gravitational potential well of those 



previously formed dark matter haloes [gOp. Dark matter haloes play a crucial role in the study of 
the formation of galaxies and the evolution of the universe, they provide the environment where 
baryons can cool down to form stellar systems, which finally are the objects we can observe in the 
universe and are the ones responsables for much of the cosmic radiation field we can account for. 

The properties of dark matter haloes can be characterized with several different parameters, 
one of the most important ones is the concentration parameter which accounts for the shape of the 
mass density profile and therefore has important implications in determining the properties of the 



galaxies forming inside [12], as well as it has direct influence in the flux of radiation produced 



during the annihilation process of dark matter in the center of the halo [0]. Previous works have 



been addressed in the study of the structure of dark matter haloes. [ |15[ ] (NFW) proposed that the 
characteristic density of dark matter haloes was directly proportional to the density of the universe 
at time of formation, making possible to connect today properties of the dark matter density profile 
to the halo formation history and to the evolution of the expanding universe. [||] and [§] have stud- 
ied the mass dependence of the concentration parameter, they found a power law mass dependence 
of the concentration parameter that scales in time with scale factor. [|19|| found a strong correlation 
between the concentration parameter and the mass accretion history of haloes, and confirmed the 
claims of NFW and ^ where the concentration of the halo is related to their definition of time of 
formation of the halo. investigated the power spectrum dependence of the concentration pa- 
rameter, re-addressed the problem of the properties of the halo mass distribution and 
studied the connection between the mass accretion history (MAH) and the concentration parameter 
using a large suite of simulations. They found a relation between the scale length of the halo and 
the mass interior to that radius M^. Using that relation in his model for the MAH and an appro- 
priated choice of the time of transition between the two different modes of accretion enabled them 
to model the mass and redshift dependence of the concentration parameter. [^] and [11] have 



studied the cosmology and mass dependence of the concentration, shape and spin parameters at 
z=0 and revised the models of NFW and [^]. In [13 1 we studied the mass and redshift dependence 



of spin, shape and concentration parameters for a WMAP5 cosmology and found the evolution of 
the concentration parameter to be due to the evolution of the inner halo mass distribution (inside 
Tg) and suggest that this evolution can be modeled as a spherical perturbation growing in the inner 
region of the halo. This description gives physical insight in to the understanding of the mass and 
redshift dependence of the concentration parameter. 



In this work we present in more detail those results presented in [ 13] and focus on the redshift 



dependence of the concentration parameter. We show our advances in the implementation of those 
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ideas in the construction of a model able to make predictions for the values of the concentration pa- 
rameter for different masses and redshifts, this may be important for the modeling of the properties 
of galaxies as well as for the study of sources of annihilation of dark matter in the early universe. 

2. Simulations and methods 
2.1 Simulations 

All simulations in this work have been performed with PKDGRAV, a tree code written by 



Joachim Stadel and Thomas Quinn [17]. The code uses spline kernel softening, for which the 
forces become completely Newtonian at 2 softening lengths. Individual time steps for each particle 
are chosen proportional to the square root of the softening length, e, over the acceleration, a: 
At I = Tj ^ ejai. Throughout, we set T] = 0.2, and we keep the value of the softening length constant 
in comoving coordinates during each run. The physical values of e at z = are listed in Table 
[T[ Forces are computed using terms up to hexadecapole order and a node-opening angle Q which 
we change from 0.55 initially to 0.7 at z = 2. This allows a higher force accuracy when the mass 
distribution is nearly smooth and the relative force errors can be large. The initial conditions are 
generated with the GRAFIC2 package The starting redshifts Zi are set to the time when the 

standard deviation of the smallest density fluctuations resolved within the simulation box reaches 
0.2 (the smallest scale resolved within the initial conditions is defined as twice the intra-particle 
distance). 

We have set the cosmological parameters according to the fifth-year results of the Wilkinson 
Microwave Anisotropy Probe mission WMAP5 [|], namely, a„, = 0.258, Q-l = 0.742, n = 0.963, 
h = 0.72, and as = 0.796, where n„, and Q.i are the values of the density parameters at z=0. Table |] 
lists all of the simulations used in this work. We have run simulations for several different box sizes, 
which allows us to probe halo masses covering the entire range lO"' h^^ Mq < M < lO'^/j^' Mq. 
In addition, in some cases we have run multiple simulations for the same cosmology and box size, 
in order to test for the impact of cosmic variance (and to increase the final number of dark matter 
haloes). 

2.2 Halo properties: Concentrations 

In this work we identify dark matter halos using a spherical overdensity (SO) algorithm and 
use a time varying virial density contrast determined using the fitting formula presented in [Q]. We 
include in the halo catalogue all the haloes with more than 500 particles inside the virial radius 
(A^vir > 500). 

To compute the concentration of a halo we first determine its density profile. The halo centre 
is defined as the location of the most bound halo particle (we define the most bound particle as 
the particle with the lowest potential energy, no care about binding energy is taken here), and we 
compute the density (p,) in 50 spherical shells, spaced equally in logarithmic radius. Errors on the 
density are computed from the Poisson noise due to the finite number of particles in each mass 
shell. The resulting density profile is fitted with a NFW profile: 

£M = ^ (21) 
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Name 


Box Size 


N 


nip 


£ 


Nrnin > 500 

z=0,2 


B20 


14.4 


250^ 


1.37e7 


0.43 


974, 1006 


B30 


21.6 


300^ 


2.68e7 


0.64 


1515, 1399 


B40 


28.8 


250^ 


1.10e8 


0.85 


1119, 993 


B90 


64.8 


6003 


9.04e7 


0.85 


13587, 12177 


B180 


129.6 


300^ 


5.78e9 


3.83 


2300, 510 


B300 


216.0 


400^ 


1.13elO 


4.74 


5840, 707 


B3002 


216.0 


400^ 


1.13elO 


4.74 


5720, 766 



Table 1: Table of simulations used in this work. Note that the name of the simulation is related to the box 
size in units of Mpc. represents the number of total particles in the box. £ represents the force softening 
length in units of kpc /i^' and the last column gives the number of haloes with more than 500 particles at 
z = and z = 2. Masses of particles are in units of /i^ ' Mg and box sizes in units of Mpc /i^ ' , with h = 0.72. 

During the fitting procedure we treat both and 5c as free parameters. Their values, and 
associated uncertainties, are obtained via a minimization procedure using the Levenberg & 
Marquardt method. We define the r.m.s. of the fit as: 

1 ^ 

Prms = £ (In pi - Inpm)^ (2.2) 

where is the fitted NFW density distribution. Finally, we define the concentration of the halo, 
Cvir = Ryir/r^, using the virial radius obtained from the SO algorithm, and we define the error on 
logc as (a,,,/rs)/ln(10), where Or^ is the fitting uncertainty on r^. 

3. Results 

In figure |l| we show the median Cyir — Myh- relation for relaxed haloes in our sample at different 
redshifts. Haloes have been binned in mass bins of 0.4 dex width, the median concentration in each 
bin has been computed taking into account the error associated to the concentration value (see ??, 
and M08). In our mass range the Cyir — Myir relation is well fitted by a single power law at almost 
all redshifts. Only for z = 2 we see an indication that the linearity of the relation in log space seems 
to break, in agreement with recent findings by [^. 

The best fitting power law can be written as: 

log(c) =fl(z)log(Mvir/[/j-^M0])+b(z) (3.1) 

The fitting parameters a{z) and b{z) are functions of redshift, the evolution of a and b can 
be itself fitted with two simple formulas that allow to reconstruct the Cvk — Myk relation at any 
redshifts: 

a{z) =wz-m (3.2) 
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Figure 1: Mass and redshift dependence of the concentration parameter The points show the median of the 
concentration as computed from the simulations, averaged for each mass bin. Lines show their respective 
linear fitting to eq. 3.1. 
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a 



■ + ■ 



(z+7) {z+rY 



(3.3) 



Where the additional fitting parameters have been set equal to: w = 0.029, m = 0.097, a = 
-110.001, j8 = 2469.720 and y = 16.885. This double fitting formulas are able to recover the 
original values of the halo concentration with a precision of 5%, for the whole range of masses 
and redshifts inspected. It has been shown by that using Nvir between 100 and 400 particles is 
enough to get good estimates for the properties of halos, nevertheless in order to look for system- 
atics we re-computed Cvir varying the minimum number of particle inside /?vir> using 200, 500 and 
1000 particles. No appreciable differences (less than 2%) were found in our results for the median. 

As can be seen from our results, the mass and redshift dependence of the concentration param- 
eter is considerably different from a simple scaling with redshift |^]. As has been already noted by 
[|l9|] and [21],p3|] the evolution of the concentration parameter is strongly coupled to the growth 
history of the halo, we will follow this idea to study the physical mechanism behind its evolution. 
To do so, we built merger trees for all haloes at z=0 in our simulations. This time we used haloes 
with N^ii- > 200 to be able to follow the evolution of haloes up to earlier times. As was already 
mentioned, using 200 particles inside the virial radius still gives good results in the estimation of 
the properties of the halo we are interested in and for simplicity we will only present results from 
our box B90. 

Our methodology is simple, since we define the concentration parameter Cvk as the ratio be- 
tween the virial radius /?vir and the scale length radius rg, we study the evolution of Cyir tracking 
the evolution of these quantities along the merger tree of haloes. It is reasonable to think that the 
properties of halos we observe at z=0 (and at any redshift) are inherited from the most massive 
progenitor, so we only trace back the properties of the halo along this branch of the tree. 

First we explore the relation between Mg, defined as the mass inside rg, and r^. Figure shows 
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Figure 2: (Left) Relation between Ms and r^ for haloes at different redshift. For all redshifts the mean slope 
equals 0.38 with very low scatter, while the normalization is clearly redshift dependent. The scatter in the 
relation is mostly due to the intrinsic scatter in the values of r^. (Right) Relation between Myir and Ms. 
Clearly both quantities are directly related with a weak redshift dependence in the normalization. 



clearly the relation between both quantities, if Ps(z) is the mean mass density inside at a given 
redshift z, then 

,3.4, 



if Ps(z) where set to be mass independent one could verify that the power index a in |3 
equals exactly 1/3, any deviation from that value would be due to a mass dependence on Ps(z)- 
The analysis on our data shows robustly that a ^ 0.38 close but clearly different from 1/3, which 
may imply that Ps{z) has a weak dependence on Mg. This behavior can be compared with the 
one expected from the relation between /?vir and Myir where the power index of the equivalent 



relation 3.4 is exactly 1/3 without any mass dependence on the mean density of the halo, which 
is constrained by the values of the overdensity contrast predicted by spherical collapse and the 
critical density of the universe. Finally, Figure || also shows that Ms and Myir are related, so, any 
dependence of Ps (z) with Ms will translate in to a dependence with Myir. 

This similarity in the relations between with Ms and /?vir with Myir suggest that one may 
think in an analogous treatment for rg and Ms as is done for /?vir and Myir- Such an analogous 
analysis is supported by the results of the time evolution of the radial scale length rs and virial 
radius /?vir along the merger trees. Figure |] shows the time (redshift) evolution of /?vir and r^. As 
it can be seen in the figure, the behavior of both quantities show a similar trend. They grow with 
decreasing redshift, reach a maximum and then start to decrease. The time at which the maximum 
is reached depends on the final virial mass of the halo, being the low mass ones the firsts reaching 
that point. 

The similarities observed in figures ^ and || suggest that one can describe the evolution of the 
inner region of the halo (the one enclosed by Tj) in a similar way as done for the outer halo. We 
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Figure 3: (Left-Top) Time evolution of the virial radius as computed from the averaged merger histories for 
haloes with final mass of 1.68 X 10' ^ 1.0 X 10'2 and 6.1 x IO'^/j-IMq. (Left-Bottom) Time evolution of 
the averaged scale length rs of the dark matter haloes in the same mass bins. (Right) As (z) as a function of 
redshift for three different mass bins computed as the ratio Ps (z) / Pvir (z) along the merger tree. 



assume the inner region of the halo to be a perturbation of mean density Ps (z) that evolves within the 
background of mean density Pviv(z) = Av,>(z)pc (z). In analogy with the spherical collapse model 
we want to look for the evolution of the density contrast of this perturbation: (z) = Ps (z) / Pvir {z) ■ 
This inner density contrast is well described by the following formula: 



^s(z) = (3-5) 
z + e{M) 

where A = 50 and e{M) = 0.3975 log (Mvir(z = 0)/[/j"^ M©]) -4.312 best reproduce our data. 



Equation implies that: i) Ps > Pvir at all redshifts, ii) Ag is a growing function of the redshift, 
implying a fast growth of the inner density with respect to the mean density of the halo and iii) 
As depends on the final mass of the halo, and it has lower values for high mass haloes. This mass 
dependence of As is also justified from the analysis of figure ^ where it was shown that the mean 
density of the inner region of the halo must be mass dependent. 

One can use all of the previous results to try to reproduce the mean values of the concentration 
parameter as a function of time for different halo masses. Figure H shows the result of modeling 



Cvir after modeling /?vii and rs with the use of eq For that Mvir(z) was modeled as a damped 



exponential law ([0], [|T3]]) and Ms was assumed to follow a power law with Myir with appropriated 
values for the power law index and normalization. 

Although the model works quite well in describing the time dependence of the concentration 
parameter, it is clear that it works much better for the high mass than for the low mass regime, 
although for the low mass haloes the model fits very well at high redshift. It is clear that the 
model performs better for halos in regimes where the nonlinear effects still are not so strong to be 
predictable by the way spherical collapse approach is implemented in this work. 
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Figure 4: Data and model of the redshift dependence of the concentration parameter Data points are 
computed as the ratio between the mean values of /?vir and rs in the respective mass bins along the branch of 
the most massive progenitors of the tree, the solid lines are the prediction of the model. 

4. Conclusions 

We present results of the study of the mass and redshift dependence of the concentration pa- 
rameter of dark matter haloes. In our mass and redshift range the Cyir — Myir relation always follows 
a power law behavior. We confirmed that the redshift dependence of such relation is more complex 
than a simple (1 + z)^' scaling as proposed by with both the normalization and the slope of the 
relation changing with cosmic time. We also found that for increasing redshifts (z 2) the power 
law behavior seems to break, in agreement with recent studies (e.g. [^). Thanks to our multiple 
box simulations we tested our results against resolution effects and find them to be stable once a 
sufficient large number of particles is used A^vir > 500. 

In order to improve our understanding on the redshift evolution of the Cyir — Myir relation 
we look at the individual evolution with time of and /?vir- Both these length scales grow with 
decreasing redshift until a maximum is reached, then they start to decrease towards z = 0. There is a 
clear analogy between the collapse of a hnear perturbation and the behavior of and /?vir- We found 
that we can model the evolution of the inner part of the halo as a decoupled spherical perturbation 
growing inside the central region of the halo. The temporal offset between the "turning points" 
of the perturbations associated with rg and /?vir is able to explain the observed redshift evolution 
of the Cvir — Mvir relation. Using this model we presented first results showing that the model we 
propose may be used to predict the evolution of the halo mass density profile for arbitrary mass and 
redshifts, which will have important applications in the modeling of the mass distribution of halos 
in redshift or mass regimes where simulations have limited resolution. 
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